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The excitation of water in the S140 photon dominated region. 
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Abstract. We consider the excitation of water in the Photon Dominated Region (PDR) S140. With the use of a three- 
dimensional escape probability method we compute the level populations of ortho- and para-H 2 up to ~ 350 K (i.e., 8 levels), 
as well as line intensities for various transitions. Homogeneous and inhomogeneous models are presented with densities of 10 4 - 
10 5 cm -3 and the differences between the resulting intensities are displayed. Density, temperature, and abundance distributions 
inside the cloud are computed with the use of a self-consistent physi-chemical (in)homogeneous model in order to reproduce 
the line intensities observed with SWAS, and to make predictions for various lines that HIFI will probe in the future. 
Line intensities vary from ~ 10~ 13 erg cm -2 s _1 sr _1 to a few times 10~ 6 erg cm -2 s _1 sr _1 . We can reproduce the intensity for 
the 1 io — > loi line observed by the SWAS satellite. It is found that the 2 i2 — > loi line is the strongest, whereas the 3i 2 — > 2 2 i line 
is the weakest, in all the models. It is found that the li — > loi line probes the total column, while higher excitation lines probe 
the higher density gas (e.g., clumps). 
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1. Introduction 

The study of water is of great interest (e.g., Ceccarelli et al. 
1996). Not only because H2O possesses a very large number of 
strong far-infrared (FIR) transitions, but also because it plays 
a crucial role in the chemistry of molecular clouds (Bergin et 
al. 1998) and can act as an important coolant in dense molecu- 
lar clouds (Neufeld et al. 1995), in shocks (Draine et al. 1983; 
Neufeld & Melnick 1987) and circumstellar outflows (Chen & 
Neufeld 1995). Because of its high abundance and easy excita- 
tion in warm interstellar and circumstellar environments, water 
is a powerful tool to probe astrophysical conditions in a broad 
variety of sources, from protostars to molecular clouds. 
Nonetheless, the large attenuation by terrestrial H2O molecules 
has made ground-based observations nearly impossible. 
Therefore, observations from space-born facilities are needed 
(e.g., Neufeld et al. 1999; Truong-Bach et al. 1999; Wright et 
al. 2000). 

The launch of the Infrared Space Observatory (ISO) has made 
it possible to look for emission and absorption lines in warm (~ 
100 K) molecular gas near young stellar objects (van Dishoeck 
et al. 1999), while the Submillimeter Wave Astronomy Satellite 
(SWAS) has probed H2O and O2 residing in cold molecu- 
lar clouds (Melnick et al. 2000), as has Odin (Hjalmarson et 
al. 2003). The future launch of Herschel, with on board the 
Heterodyne Instrument for the Far Infrared (HIFI), will pro- 
vide even better information on the physical and chemical con- 



ditions in molecular clouds by means of observing many water 
lines with higher angular resolution and sensitivity. 
In this paper we discuss the excitation of water in PDRs. As a 
testcase we take SI 40 because it has been probed extensively 
in the past with SWAS and will be examined in the future with 
HIFI. The results presented here are obtained with the use of a 
newly developed escape probability method described in Sect. 
2 combined with the inhomogeneous numerical code of Spaans 
(1996), expanded in Spaans & van Dishoeck (1997). Section 3 
describes briefly S140. Section 4 is used to describe the results 
from the 3D- (in)homogeneous models. We summarize and dis- 
cuss the paper in Sect. 5. 



2. The escape probability method 

We calculate the transfer of line radiation of ortho- and para- 
H2O in a (in)homogeneous 3 -dimensional spherical cloud by 
use of an escape probability approximation. Let ni(x, y, z) cm -3 
be the population density of the z th level in point (x, y,z). The 
equations of statistical equilibrium can then be written as: 



/ 

ni(x,y,z)^Rij(x,y,z) 



^\ij(x,y, z)Rji(x,y, z), 



(1) 



Send offprint requests to: D.R.Poelman@astro.rug.nl 



where / is the total number of levels included. Rij(x, y, z) is ex- 
pressible in terms of the Einstein A t j and B t j coefficients, and 
the collisional excitation (i < j) and de-excitation rates (i > j) 
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Rij(x,y,z) 



\Aij + Bij < Jij > + Cij(x, y, z), (i > j) 



\Bij < > + Cij(x,y,z), 



(i<j) 



(2) 



<Jij(x, y, z)> is the mean integrated radiation field at frequency 
Vij corresponding to the transition from level / to j at position 
(jc, y, z) in the cloud, given by: 

<Jij(x,y,z)> = [1 - Pij(x,y,z)]Sij(x 9 y 9 z) 

+Pi j (x,y,z)B ij (v ij ,T BB ), (3) 

where ^ ; (x,y,z) is the probability for a photon to escape and 
S//(x,y,z) is the source function, to be addressed later. The 
background radiation Bij(vij,T BB ) consists of two terms: the 2.7 
K microwave background and the infrared emission of dust at 
a temperature and with an emission optical depth Td(v^): 

Bij(y ij9 T BB ) = %J = 2.7K) 

+r d (y i j)B(y ij ,T d ), (4) 

where r^ivij) = Tioo^mQOO fim/A). We adopt a value for Tiooyum 
of 0.001 (Hollenbach et al. 1991). If all the photons escape 
(e.g., Pi] = 1), then <Jij(x,y,z)> is ^(v/^Tbb); if none es- 
cape (e.g., Pij = 0) <Jij(x,y,z)> is the local source function 

Sij(x,y, z). 

The source function reads as 



Sij(x,y,z) = 



2hv 3 . 



rij(x,y,z)gi 



1 



(5) 



_ni(x,y,z)gj 

where gi and gj are the statistical weights of levels / and j, re- 
spectively. 

Since the set of / statistical equilibrium equations is not inde- 
pendent, one equation has to be replaced by the conservation 
equation 

/ 

nx = ^Jij, (6) 

7=0 

where n x is the number density of species x in all levels. 
In a sphere, the probability that a photon emitted in the transi- 
tion from level i to level j at position (x, y, z) along a direction 
k escapes the cloud is given by 

1 - Qxp(-Tij(x,y,z,k)) 



Pij(x,y,z, k) = 



(7) 



Tij(x,y,z, k) 

with nj the optical depth in the line. Note that the probability 
for a photon to escape in a point, p(x,y,z), equals the sum of 
the escape probabilities over all directions, e.g., fiij(x, y, z) = 
J3ij(x,y,z,k). The number of directions is arbitrary, but a 6-ray 
approximation is implemented for the 3D models. Note that for 
a photon emitted in a plane-parallel medium, the probability to 
escape is given by the expression 

fN 1 - exp(-3Tij(x,y 9 z,k)) 
Pij(x,y,z,k) = - 



3Tij(x 9 y,z 9 k) 



(8) 



The optical depth averaged over the line, in direction k over a 
distance s = s 2 (x 2 ,y2,Z2) - si(x\,yi,zi) is given by 

S-S2 



Tij(x,y,z,k) = —r — 1 

8ttv?. J Av d [mgf 



ds, 



(9) 



with Avd = (v^ h + v^ urb ) 1/2 the velocity dispersion. Turbulence 
in PDRs can be combined with the thermal speed of the gas 
into a Gaussian distribution. The resultant Doppler profile is 



(f(x) 



-x 2 

e x 



(10) 



where x = (v - v^)/Avd. The Doppler frequency width Avd is 
given by Avd = Avd v -f. The turbulent dispersion, v tur b, is larger 
than the thermal velocity of the gas, v t h, and dominates the line 
broadening in PDRs. 

Because interpretation of observed water lines requires knowl- 
edge of collisional excitation rates responsible for line forma- 
tion, it is necessary to work with the most accurate collisional 
rate data available. The collisional rates are taken from Green 
et al. (1993) for collisions involving up to 8 levels for inelastic 
collisions between H 2 and He. To account for the different 
reduced mass of the collisional partner (e.g., H2), the rate co- 
efficients are scaled by a factor of 1.348. Rate coefficients of 
para- and ortho H 2 owing to collisions with ortho- and para- 
H2 at kinetic temperatures from 20 to 140 K are taken from 
Phillips et al. (1996). In addition, rate coefficients for H2O in 
collisions with P-H2 and 0-H2 for the low temperature regime 
(5-20 K) are taken from Grosjean et al. (2003) and Dubernet 
& Grosjean (2002), respectively. Clumps and their edges dom- 
inate the total water emissivity a few magnitudes of extinction 
into the PDR. The ambient electron abundance in these regions 
does not exceed 3 x 10" 7 -10 -6 in our models. Using the H2O- 
electron collision rates presented by Faure et al. (2004) we find 
that the contribution by electrons to the total collisional excita- 
tion rate is less than 1 % and is therefore ignored in our calcu- 
lations. Note that when the electron abundance exceeds 10" 5 , 
this contribution is equal to the collisional rates of H2O with 
para-H2, and consequently has to be taken into account. 
In PDRs, the ortho-to-para ratios of H2 have typical values in 
the range 1.5-2.2 (Hasegawa et al. 1987; Ramsay et al. 1993; 
Chrysostomou et al. 1993; Hora & Latter 1996; Shupe et al. 
1998), and are not constant because of gas temperature varia- 
tions within the PDR. Therefore we adopt the expression for 
the ortho-to-para ratio (OPR), in thermal equilibrium, defined 
by 



OPR = 



(2/ + 1)2(27+ l)exp(-^f^) 
(2/ p + 1) 2(2/ + 1) exp(- g^g^ ) 



(11) 



where I and I p are the total nuclear spin, according to the hy- 
drogen nuclear spins being parallel (7 = 1, TT) or anti-parallel 
(7p = 0, TD- The sum in the numerator (denominator) extends 
over all ortho (para) levels (7, K c ) (Mumma et al. 1987). A 
similar expression holds for the ortho- and para-I^O OPR, for 
which the statistical equilibrium (OPR = 3) value is attained at 
temperatures exceeding ~ 60 K, but will differ from the statis- 
tical equilibrium value for lower temperatures. 
This escape probability method is the first step towards a fully 
3D-Monte Carlo radiative transfer code that calculates line in- 
tensities, as well as line profiles in different astrophysical envi- 
ronments. 
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3. The molecular cloud S140 

The molecular cloud lies at a distance of ~ 900 pc (Crampton 
& Fisher 1974) and is illuminated from the south-western side 
by the BOV star HD 211880. The star is located about ~ T 
(1.85 pc) from the edge of the cloud. The cloud extends over 
more than 30' (8 pc) and contains a dense core in which star 
formation is occuring. 

The molecular cloud consists of two parts: (i) the extended 
molecular cloud (EMC), and (ii) the interface region near the 
south-west edge of the cloud (PDR). Here focus is on the in- 
terface region, as this region is a far more interesting object to 
probe the behaviour of water with its high H2O abundance rel- 
ative to the extended molecular cloud. 

In the past, models were neither able to reproduce the obser- 
vations of extended emission of [Ci] (Keene et al. 1985) and 
[Cn] (Stutzki et al. 1988) deep into the cloud, nor the intense 
13 CO 6-5 emission in some regions (Graf et al. 1990). This in- 
formation, together with the CO maps by Falgarone & Perault 
(1988), led to the suggestion that interstellar clouds have a very 
inhomogeneous structure, so that the ultraviolet radiation can 
penetrate much deeper into the cloud. 

Most previous models of S140 (e.g., Ashby et al. 2000) have 
assumed a spherically symmetric cloud with temperature and 
density power laws. Even though they can reproduce well the 
observed value of the intensity of the 557 GHz lio — > loi 
ground-state transition of ortho-H 2 16 0, being acquainted with 
the inhomogeneous structure of molecular clouds from obser- 
vations and in order to explain and predict more accurately var- 
ious line transitions, we make use of density, temperature and 
abundance (relative to H2) distributions as calculated with the 
inhomogeneous code of Spaans (1996). This code is described 
first in Spaans (1996) and modified in Spaans & van Dishoeck 
(1997). This is the first time an attempt is made to model the 
water emission in SI 40 in an inhomogeneous way. 
The Monte Carlo code calculates the chemical structure and 
thermal balance in PDRs, as well as the distribution of the 
water abundance and molecular hydrogen density simultane- 
ously. A chemical network including 291 reactions between 
51 species consisting of the elements hydrogen, carbon ( 12 C), 
oxygen, iron, magnesium, and PAHs (with charges -1, 0, +1) 
is used. The gas abundances for S140, are taken to be 2.0 x 
10- 4 , 5.0 x 10- 4 , 2.5 x 10- 7 and 1.3 x 10" 6 for C, O, Fe and 
Mg, respectively. For carbon this implies that about 50% of the 
solar abundance is present in the gas phase, which is compa- 
rable with the values found in diffuse clouds (Cardelli et al. 
1993). Mechanisms for converting the ultraviolet flux into gas 
heating are (1) photoelectric emission from small dust grains 
and large molecules like Polycyclic Aromatic Hydrocarbons 
(PAHs) (Bakes & Tielens 1994; Spaans et al. 1994), (2) neutral 
carbon ionization, (3) H2 photodissociation, (4) FUV pump- 
ing followed by collisional de-excitation of vibrationally ex- 
cited H2 in the two level approximation, and (5) H2 formation 
from FUV-produced hydrogen atoms on grain surfaces. In ad- 
dition, (6) cosmic-ray heating may also contribute deep into 
the cloud, and (7) 1 may heat the gas through the collisional 
de-excitation of the 3 Pi fine structure level excited by (dust) 
continuum radiation. 




Fig. 1. Gas (left) and dust (right) temperatures (K) when taking a ver- 
tical slice through the centre of the sphere in case of a homogeneous 
model. The star HD 21880 is located at the left of both images. 

For the cooling, the infrared fine-structure lines from atoms C, 
C + , O, Si, and Fe, if present in the chemical network, have 
been taken into account. The lowest 20 rotational levels of CO 
are calculated in statistical equilibrium and their contribution 
to the cooling is included. H2 rotational cooling is considered 
as well as cooling due to water through a cooling function 
(Neufeld & Melnick 1987). The contribution of water to the 
total cooling rate amounts to 10-20% around the clump den- 
sity of 10 5 cm" 3 and for temperatures of 40-60 K. The CO 
molecule dominates the cooling for these conditions. 
The computed temperature, density and water abundance dis- 
tribution is then used as input for the water line transfer code 
described previously in this paper. This approach leads to a 
more accurate treatment of the complex chemical and thermal 
structure of the PDR, than when a temperature and H2O abun- 
dance are simply imposed. 

4. Models 

This section consists of two parts. The first part describes the 
results from the homogeneous model, the second part describes 
the inhomogeneous models and the results thereof. 

4.1. Homogeneous 3D model 

The incident radiation field, taken to be /uv = 140 with respect 
to the Draine (1978) field, is consistent with the enhancement 
at the edge of the cloud near the BOV star HD 211880. An 
isotropic component, /uv = 1, of the average interstellar ra- 
diation field has also been taken into account. For the homo- 
geneous PDR model we adopt a generic mean total hydrogen 
density of nu = n(H) + 2ft (H2) = 2 x 10 4 cm -3 over a region 
0.5 pc in extent, or Ay ~ 20 mag, in agreement within a factor ~ 
1.5-2 with observations (Hayashi et al. 1985; Zhou et al. 1994). 
For areas where H2O is present, hydrogen is totally molecular 
(the main part of the UV radiation is already absorbed), and 
therefore nu = 2n(H 2 ). The grid size is taken to be 81 x 81 x 
81, corresponding to a resolution of 0.006 pc. 
The temperature distribution for the gas and dust is shown in 
Fig. 1 . The gas (dust) temperature is « 150 K (40 K) at the west- 
ern edge of the cloud, and drops to about 30 K (30 K) deeper 
into the cloud. The dust temperature stays far below the gas 
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Table 1. Homogeneous model 

A u \ transition E u \ A u \ average 61 peak average n cx b 

intensity intensity around 

peak intensity c 











[K] 


[s] 




[erg cm 2 s 1 


sr 1 ] 


[erg cm 2 s 1 


sr 1 ] 


[erg cm 2 s 1 


sr 1 ] 


[cm 3 ] 


538.2 


o- 


H 2 1 10 


-> loi 


26.7 


3.458 x 10" 


-3 


4.00 x 10 


-8 


1.00 x 10" 


-7 


9.20 x 10" 


-8 


3.5 x 10 7 


179.5 


o- 


H 2 2i2 


-> loi 


80.2 


5.593 x 10" 


-2 


5.66 x 10 


-7 


1.41 x 10" 


-6 


1.37 x 10" 


-7 


5.0 x 10 8 


108.1 


o- 


H 2 2 2! 


-> lio 


133.2 


2.564 x 10" 


-1 


8.26 x 10" 


10 


2.00 x 10" 


-9 


1.88 x 10" 


-9 


3.6 x 10 9 


180.5 


o- 


H 2 22i 


-> 2i 2 


79.7 


3.058 x 10" 


-2 


5.90 x 10" 


11 


1.43 x 10" 


10 


1.35 x 10" 


10 


6.6 x 10 8 


174.6 


o- 


H 2 3o3 


-> 2i 2 


82.4 


5.048 x 10" 


-2 


5.70 x 10" 


10 


1.04 x 10" 


-9 


8.98 x 10- 


10 


7.7 x 10 8 


260.0 


0- 


H 2 3i 2 


-> 2 2 i 


55.4 


2.634 x 10" 


-3 


1.86 x 10" 


12 


4.04 x 10" 


12 


3.78 x 10" 


12 


4.8 x 10 7 


273.2 


0- 


H 2 3i2 


-> 3 03 


52.7 


1.648 x 10" 


-2 


1.11 x 10" 


11 


2.40 x 10" 


11 


2.26 x 10" 


11 


2.8 x 10 8 



269.3 p-H 2 In 

303.4 p-H 2 2 02 ■ 
101.0 p-H 2 2 20 
398.6 p-H 2 2n 



>0 00 
►In 
►In 

> 2q 2 



53.5 
47.5 
142.6 
36.1 



1.842 x 10- 2 
5.835 x 10" 2 
2.607 x 10" 1 
7.062 x 10" 3 



1.76 x 10" 7 
1.78 x 10- 9 
1.93 x 10" 10 
3.34 x 10" 11 



4.94 x 10- 7 
4.42 x 10" 9 
5.40 x 10" 10 
6.40 x 10" 11 



4.54 x 10- 7 
4.18 x 10" 9 
5.00 x 10" 10 
5.64 x 10" 11 



1.6 x 10 8 
8.9 x 10 7 
3.8 x 10 9 

8.7 x 10 7 



a Is the average intensity over the whole cloud. 
b Critical density evaluated at 100 K. 

Is the intensity averaged over the HIFI FWHM beam size 
(http://www.sron.rug.nl/hifiscience/) for a beam centered on the peak 
emission. 




5e-08 le-07 le-12 2e-12 3e-12 




2e-07 4e-07 2e-ll 4e-ll 



Fig. 2. Distribution of the intensities (erg cm -2 s _1 sr _1 ) for various 
transitions of ortho- (blue/top) and para- (red/bottom) H 2 0. The top 
left panel depicts the 1 10 — > loi transition, the top right the 3 i2 — > 
2 2 i transition. Bottom left is In — > Ooo, and bottom right is 2\\ — > 2 02 . 
The star HD 21 1880 is located to the left of the Figures. 



temperature at the west side of the cloud, but becomes slightly 
larger or equal to the gas temperature farther into the cloud. It is 
worth noting that dust plays an important role in the excitation 



of water. It absorbs the interstellar FUV and visible radiation 
field and emits the absorbed energy in the far infrared, offering 
a source of excitation for water. Computed water abundances 
vary from ~ 10" 10 to ~ 10" 6 , relative to the hydrogen density 
with an average abundance of 1.4 x 10" 7 . At the west side 
of the cloud, where the star is located, the abundance is lowered 
due to photodissociation of water by the FUV radiation field: 

H 2 + hv — > H + OH, 
OH + hv -> O + H 

Deeper into the cloud, the radiation field is attenuated by dust, 
resulting in a higher abundance. 

At low temperatures H2O is formed mainly through the disso- 
ciative recombination of the H30 + ion: 

H 3 + + e" -> H 2 + H 

A value of 0.3 is chosen for the branching ratio leading to H2O 
in the dissociative recombination of HsO + with a rate coeffi- 
cient of 3.3 x 10 -7 (77300)- - 3 cm 3 s _1 . The branching ratio 
is consistent with the results of Vejby-Christensen et al. (1997) 
but higher than the value of 0.05 suggested by Williams et al. 
(1996). At temperatures above ~ 300 K, the formation of wa- 
ter is dominated by the neutral-neutral endoergic reaction se- 
quence: 

0+ H 2 ^ OH + H, 
OH + H 2 ^ H 2 + H 

(Elitzur & de Jong 1978; Elitzur & Watson 1978). The activa- 
tion energy barriers are overcome by the thermal energy of the 
warm gas. 

The emergent intensity integrated along a line of sight is 
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calculated in the following way: 



/ = (47T)" 1 f h{z')dz' , 

Jo 

where 

Myij) = ruAijhvijPiTij) 



(12) 



(13) 



with Hi the population density in the i level, A*/ the Einstein A 
coefficient, hv i; the energy difference between the levels i and 
j, Ti 7 - the optical depth averaged over the line, and p(jij) the es- 
cape probability at the optical depth nj of the line. 
Average intensities, peak intensities, as well as intensities av- 
eraged over the HIFI FWHM beam size at the transition fre- 
quency of the particular line for a beam centered on the peak 
emission are listed in Table 1 . Figure 2 shows maps of the pre- 
dicted distribution of the intensities for various transitions for 
ortho-, and para-H 2 as we would see them on the sky, con- 
volved with a beam size of 0.3" (-the angular extent that cor- 
responds to one gridpoint). 

To interpret the results one has to take into account the distribu- 
tion of the gas and dust temperature (T g and respectively), 
the abundance of water in the cloud, the density of the medium, 
as well as optical depth effects. Since the overall mean den- 
sity of hydrogen is 2 x 10 4 cm -3 , far below the critical density 
for all the transitions (of the order of 10 8 cm -3 or higher), the 
lines are subthermally excited. This means that collisional de- 
excitation is less important than spontaneous emission, which 
results in a distribution of the level populations deviating from 
LTE. 

The emission will originate from a depth where r ;j ~ 1. The 
place where r ;j ~ 1 varies for every transition. Optically thick 
lines, i.e., > 1, can be used to probe the physical condi- 
tions in the outer (t^- < 1) layers of the cloud. Lines which 
are optically thin can be used to probe the physical conditions 
throughout the cloud. A peak in intensity for a certain transi- 
tion will occur at places in the cloud where the temperature is 
high enough to excite the line and will depend on the column 
of H2O along the line of sight. Therefore at the edges of the 
cloud, where the column is rather low, the intensity drops. 
The 0-H2O 1 10 — > loi transition peaks near the east edge of the 
cloud where the temperature of the gas drops to ~ 60 K, the 
perfect regime to excite the line. Transitions with a higher en- 
ergy difference, e.g., 0-H2O 3 12 —> In , have their peak shifted 
towards the edge of the cloud where the temperature is higher. 
Note that HIFI, with its high angular resolution will be able to 
see this shift. Depending on the frequency of the transition, the 
beam will be ~ 5 times smaller than the PDR. 
The 0-H2O 2 12 — » loi transition is the strongest we found, with 
an average intensity of « 5 x 10" 7 erg cm -2 s" 1 sr" 1 . 






Fig. 3. Distribution of the intensities (erg cm -2 s _1 sr -1 ) for various 
transitions of ortho- H 2 using the PDR I model. The top left panel 
depicts the lio — > loi transition, the top right the 2 i2 — > loi transition. 
Bottom left is 2 2 i — > 1 10, and bottom right is 3i 2 — > 2 2i . The star HD 
21 1880 is located to the left of the Figures. 






4.2. Inhomogeneous 3D model with the EMC 

The inhomogeneous nature of molecular clouds is well es- 
tablished from observations of extended [Ci] and [Cn] (e.g., 
Keene et al. 1985; Plume et al. 1994) and has been modeled by 
various groups (e.g., Meixner & Tielens 1993; Spaans 1996; 



Fig. 4. Distribution of the intensities (erg cm -2 s _1 sr -1 ) for various 
transitions of para-H 2 using the PDR I model. The top left panel 
depicts the In — > o transition, the top right the 2 02 — > In transition. 
Bottom left is 2 20 — > In, and bottom right is 2n — > 2 02 . The star HD 
21 1880 is located to the left of the Figures. 
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Table 2. Model parameters 



Model 


size 


n c 




F 


Ric 


k 


A(H 2 0) fl 


Av d 




[pc] 


[cm" 3 ] 


[cm" 3 ] 


[%] 




[pc] 




[km s" 1 ] 


PDRI 


0.5 


4x 10 5 


1 x 10 4 


2.5 


40 


0.03 


1.1 x 10~ 8 


1.2 


PDR II 


0.5 


2x 10 5 


1 x 10 4 


5 


20 


0.03 


1.7 x 10- 8 


1.2 


PDR III 


0.5 


6x 10 4 


6x 10 3 


20 


10 


0.03 


2.1 x 10" 8 


1.2 


EMC 


10 


3000 


300 


25 


10 


0.6 


6.6 x 10" 10 


0.7 



a Volume averaged water abundance. 



Table 3. Inhomogeneous PDR models 



Model 


transition 


average 


peak 


average 








intensity 


intensity 


around 












peak intensity 61 








[erg cm -2 s _1 sr _1 ] 


[erg cm -2 s _1 sr _1 ] 


[erg cm -2 s _1 sr _1 ] 


PDRI 




^ 1()1 


1.06 x 10" 7 


1.92 x 10" 6 


2.20 x 10- 7 




u ri2 w z*YL 




6.16 x 10~ 7 


6.56 x 10~ 6 


1.39 x 10~ 6 




u ni w z<2i 


^ J-IO 


5 65 x 1 0" 9 


1.71 x 10" 7 


9.12 x 10" 9 




o-H 2 22i 


-> 2i2 


4.05 x 10" 10 


1.23 x 10" 8 


6.51 x 10" 10 




O-H2O 3o3 


^2i2 


4.67 x 10- 9 


1.20 x 10- 7 


8.00 x 10- 9 




0-H 2 3i2 


^2 2! 


1.31 x 10- 11 


5.38 x 10- 10 


2.81 x 10- 11 




0-H 2 3i2 


— * -^03 


7.84 x 10- 11 


3.20 x 10- 9 


1.68 x 10- 10 




p-H 2 o 111 


— > Ooo 


2.50 x 10" 7 


3.69 x 10" 6 


7.16 x 10- 7 




P-H2O 2q2 


-> 111 


2.70 x 10- 8 


8.08 x 10- 7 


5.83 x 10- 8 




p-H 2 2 20 


-> 111 


8.80 x 10- 10 


2.29 x 10- 8 


1.72 x 10- 9 




P-H 2 2n 


^2 02 


2.69 x 10- 10 


7.80 x 10- 9 


4.39 x 10- 10 


PDR II 


0-ri2^ lio 


~~ * loi 


7.25 x 10" 8 


6.72 x 10" 7 


1.55 x 10- 7 




n N,D 9,o 
v' rn w Z42 


^ J-01 


7.30 x 10~ 7 


5.71 x 10~ 6 


1.81 x 10" 6 




n H^n 9^ 1 
U-JI2U z<2l 


v 1 ^ 

^ J-IO 


1.79 x 10" 9 


2.59 x 10" 8 


3 09 x 1 0~ 9 




o-H 2 2 2 i 


-> 2l9 


1.27 x 10" 10 


1 85 x 10" 9 


2.18 x 10" 10 




o-H 2 3 03 


^2i2 


1.48 x 10~ 9 


2.24 X 10~ 8 


2.61 x 10~ 9 




0-H 2 3i2 


^2 2 i 


4.28 x 10- 12 


7.84 X 10- 9 


1.06 x 10- 11 




0-H 2 3i2 


3 3 


2.54 x 10- 11 


4.65 X 10- 10 


6.34 x 10- 11 




P-H2O In 


^Ooo 


2.52 x 10- 7 


2.47 x 10" b 


8.64 x 10- 7 




p-H 2 2 02 


In 


7.72 x 10- 9 


1.50 x 10- 7 


1.72 x 10- 7 




p-H 2 2 20 


In 


3.53 x 10- 10 


5.93 x 10- 9 


7.84 x 10- 10 




P-H 2 2n 


-> 2 02 


8.35 x 10- 11 


1.51 x 10- 9 


1.81 x 10- 10 


PDR III 


0-H2O lio 


-> loi 


2.66 x 10" 8 


2.17 x 10" 7 


1.07 x 10- 7 




0-H 2 2i2 


->loi 


3.94 x 10- 7 


2.95 x 10- 6 


1.58 x 10- 6 




o-H 2 2 2 i 


lio 


1.97 x 10- 10 


1.60 x 10- 9 


5.93 x 10- 10 




o-H 2 2 2! 


-> 2l2 


1.41 x 10- 11 


1.14 x 10- 10 


4.23 x 10- 11 




o-H 2 3 03 


^2i2 


8.81 x 10- 11 


3.82 x 10- 10 


1.04 x 10- 10 




0-H 2 3i2 


^2 2 i 


2.06 x 10- 13 


1.54 x 10- 12 


5.31 x 10- 13 




0-H 2 3i2 


3 3 


1.23 x 10- 12 


9.12 x 10- 12 


2.68 x 10- 12 




P-H2O In 


-> Ooo 


1.48 x 10" 7 


1.41 x 10" b 


6.88 x 10- 7 




P-H 2 2 02 


111 


4.55 x 10- 10 


3.95 x 10- 9 


1.08 x 10- 9 




p-H 2 2 20 


111 


6.93 x 10- 11 


9.51 x 10- 10 


2.75 x 10- 10 




P-H 2 2n 


-> 2 02 


4.81 x 10- 12 


2.72 x 10- 11 


7.39 x 10- 12 



a Is the intensity averaged over the HIFI FWHM beam size 
(http://www.sron.rug.nl/hifiscience/) for a beam centered on the peak 
emission. 



Stoerzer et al. 1996). It is therefore natural to explore the ex- 
citation of H 2 in models incorporating high density clumps 
embedded in a low density interclump medium. 
In the case of an inhomogeneous distribution, the radiative 
transfer is solved in three dimensions ona81x81x81 grid. 
In this section we not only cover the excitation of water in the 



PDR but also in the EMC, in order to fully interpret the results. 
We define the molecular cloud core as that part of the EMC 
that encloses the PDR and lies in one SWAS beam. What is left 
outside the molecular cloud core is denoted as the EMC. 
The clumpy models have three free parameters: (i) the volume 
filling factor F, which gives the fraction of the total volume 
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depth in the cloud [pc] 



Fig. 5. Change in intensity along a W-E cut for the ortho li — > loi 
(blue/top), and ortho 3\ 2 — > 2 2 i (red/bottom) transition as seen in Fig. 
3. The 3 12 — > 2 2 i transition is multiplied by 1000. The intensity, for the 
1 io — > loi rises until it reaches its maximum near the east side of the 
PDR and slightly levels off towards the edge. The 3 i2 — > 2 2 i transition 
reaches its peak near the west side of the PDR and levels off towards 
the east side. 



Table 4. Inhomogeneous EMC model 



Model 


transition 


average 






intensity 






[erg cm -2 s _1 sr _1 ] 


EMC core 


o-H 2 lio -> loi 


1.37 x IO" 11 




o-H 2 2i2 -> l i 


5.45 x IO" 11 




o-H 2 2 2i -> lio 


2.78 x 10~ 14 




o-H 2 2 2 i -> 2i2 


1.98 x IO" 15 




o-H 2 3 03 -> 2i2 


2.90 x IO" 14 




o-H 2 3i2 -> 2 2i 


1.53 x IO" 16 




o-H 2 3i2 -> 3 03 


9.17 x IO" 16 




p-H 2 In -> 00 


2.49 x IO" 11 




p-H 2 2q2 ^ In 


6.04 x IO" 14 




p-H 2 2 20 -> In 


5.19 x IO" 15 




p-H 2 2n -> 2 02 


2.23 x IO" 15 



Table 5. Antenna temperatures of the o-H 2 lio — > loi transition 



Model 


antenna temperature 

n 

[K] 


Homogeen 


0.08 


PDR I 


0.25 


PDR II 


0.17 


PDR III 


0.07 


EMC core 


3.6 x IO" 4 


EMC 


3 x 10- 4 -3 x IO" 3 



that is occupied by the clumps, (ii) a clump size / c , which fixes 
the extinction through an individual clump, and (iii) the clump- 
interclump density ratio R lc . For the PDR a fixed incident radi- 
ation field has been taken Ijjy = 140 with respect to the Draine 
(1978) field, as well as a mean total hydrogen density <nu> = 2 
x 10 4 cm -3 , as in the case of the homogeneous PDR. The mean 
total hydrogen density <nu> obeys <nu> - Fn c + (1 - F)n lc . 
Here n c and n lc indicate the clump and interclump density, re- 



spectively. The enhancement factor for the EMC decreases to 
Ijjy ~ 30-50 because of geometric dilution. The mean total hy- 
drogen density <nu> is taken to be 10 3 cm -3 over a region 
10 pc in extent. The average hydrogen density of the EMC is 
not well constrained by observations, but the adopted value ap- 
proximates those derived from the total 13 CO column density 
divided by the size of the cloud (Blair et al. 1978; Plume et al. 
1994). 

Table 2 gives the parameters used for the different models. 
Three models have been run for the PDR and one for the EMC. 
The PDR clump densities range from 3xl0 4 cm -3 to 2xl0 5 
cm -3 , whereas interclump densities vary between 3xl0 3 cm -3 
and 5xl0 3 cm -3 . Different volume filling factors were incor- 
porated, ranging from 2.5 % to 20 %. In the PDR models, the 
clump size has been taken constant, / c = 0.03 pc. The clump 
size in the EMC is taken to be 20 times larger than the clump 
size in the PDR, equal in proportion to the ratio in scale size. 
The velocity dispersion is taken to be 1.2 km s" 1 in the PDR 
(Zhou et al. 1994), and 0.7 km s" 1 in the EMC, assuming a 
weaker turbulent velocity component in comparison with the 
Vturb of the PDR. The larger PDR velocity width simply reflects 
our assumption that the embedded sources inject some addi- 
tional turbulent motions (IRS 1-3, Blair et al. 1978; Beichman 
et al. 1979; Evans et al. 1989). This seems reasonable given that 
star-forming molecular cores tend to have larger (non-thermal) 
velocity widths than quiescent cores. 

The clumps are distributed randomly in the interclump 
medium. The clump and interclump gas join smoothly with a 
density gradient. Due to the dumpiness the radiative transfer 
will differ from that of the homogeneous model. The UV radi- 
ation penetrates deeper into a clumpy region than into a cloud 
in which the gas is homogeneously distributed, leading to a 
lower water abundance because of enhanced photodissociation 
(Spaans & van Dishoeck 2001). The average water abundance 
in comparison with the homogeneous model drops by a fac- 
tor ~ 10. However, the resulting intensities depend sensitively 
on the optical depth, density and temperature structure along 
the line of sight. The inhomogeneity also affects the tempera- 
ture structure, resulting in enhanced temperatures deeper into 
the interclump medium of the cloud, and lower clump temper- 
atures. The result is an increase in the volume averaged gas 
temperature (~ 78 K for PDR I), compared with that of the ho- 
mogeneous PDR model (~ 52 K). However, the mass-weighted 
average gas temperatures in the inhomogeneous models are 88 
K, 70 K and 46 K for PDR I, PDR II and PDR III, respectively, 
and 52 K for the homogeneous model. This means that most of 
the mass of the cloud is located in cold clumps. 
Intensities are computed in the same way as in the homoge- 
neous model using Eq. (12) and (13). We calculate the total 
intensity along each line of sight. Average and peak intensi- 
ties for various transitions for the PDR are shown in Table 3, 
and for the molecular cloud core, i.e., that part of the EMC ly- 
ing at the southwest of the cloud which encloses the PDR, in 
Table 4. Fig. 3 and 4 show maps of the predicted distribution 
of the intensities for various transitions, and Fig. 5 shows the 
change in intensity along a W-E cut for two ortho- H2O transi- 
tions. Note that the vast majority of the EMC is not contained 
in the SWAS beam. Snell et al. (2000) mapped a larger region 
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of SI 40 with SWAS, including several positions on the molec- 
ular cloud. They find no significant 557 GHz water emission 
residing in these areas (e.g., outside the molecular cloud core), 
in agreement with our results. The contribution of the EMC 
(core) to the total signal is negligible in comparison with the 
PDR. The minor contribution of the EMC (core) is due to the 
low gas density, which results in little excitation of the H2O 
molecule. Freeze-out of water is not required to explain the 
SWAS non-detections away from the PDR. However, the de- 
rived physical conditions indicate that freeze-out of water does 
occur. From Bergin et al. (1995) we find a depletion factor of 
about 2-3 (see Section 5). 

The intensity for the 0-H2O 1 10 — > loi line in the homogeneous 
PDR model is a factor of a few lower than in PDR I and PDR II, 
and equals that of PDR III. The strongest transition we find is 
the o-H 2 2 12 — > loi in all models. This transition has an upper 
state energy of ~ 80 K above ground level. The gas tempera- 
tures in the models are sufficient to excite this transition. All 
other transitions are a few orders of magnitude weaker than the 
two lowest transitions. To excite these transitions much higher 
temperatures are needed. For para-H20, the lowest transition 
In — > Ooo is found to be the strongest in all the models. 
From the average intensities, antenna temperatures are derived 
for the 557 GHz transition as seen in Table 5, corrected for 
beam dilution in the case of the PDR with a beam dilution fac- 
tor of 2. Note that the correct beam dilution factor in case of 
the PDR is uncertain. Since the antenna temperatures scale in- 
versely in proportion to the beam dilution factor the results pre- 
sented in Table 5 will differ when another beam dilution factor 
is incorporated. 

The total antenna temperature is the weighted sum of the PDR 
and EMC contributions. PDR I and PDR II in Table 3 yield, 
when combined with Table 4, a value consistent with the ob- 
servation of SWAS, i.e., « 0.25 K, (Ashby et al. 2000; Snell 
et al. 2000) for the lowest lying transition of ortho-H 2 0. At 
face value PDR I is in closest agreement with the observations. 
Maps for the various transitions (line, not line + continuum) of 
the PDR I model are plotted in Figs. 3 and 4. It is seen that the 
0-H2O 1 10 — > loi ground state transition peaks near the eastern 
edge of the cloud. The bulk of the emission resides in the indi- 
vidual clumps. These clumps are low-temperature- water reser- 
voirs. At the west side of the cloud the water abundance drops 
due to dissociation, still, for higher transitions, we see that the 
majority of the emission is shifted towards the western edge of 
the cloud. The H2O emission arises mainly from the warm, ~ 
60-110 K, clump edges, where the gas temperature is higher 
than in the interiors of the clumps. Conversely, the clump edge 
water abundance is a factor of ~ 10 lower compared to the in- 
sides of the clumps, but still a factor of ~ 10 higher than in 
the interclump medium. The 557 GHz line traces mostly the 
total column of water, while it is evident from Table 3 that the 
higher excitation lines differ from one model to the other more 
strongly. 

Again, HIFI will be able to confirm the predicted shift towards 
the western edge of the higher water transitions because of its 
good angular resolution. Also, HIFI can constrain the inhomo- 
geneity of the medium through water line observations. 



5. Summary and Discussion 

In this paper we have examined the influence of inhomogeneity 
on the excitation of water in molecular clouds. Inhomogeneity 
causes photons to penetrate much deeper into the cloud lead- 
ing to enhanced photodissociation of water. This gives rise to a 
lower water abundance. However, calculated intensities of PDR 
I and PDR II are a factor a few higher for the case of a ho- 
mogeneous distribution. PDR III is in closest agreement with 
the intensities of the homogeneous model. This is not surpris- 
ing since the parameters of PDR III resemble the homogeneous 
model most. This leads to the conclusion that the major part of 
the water line intensity emanates from the warm high density 
clump edges. 

The total antenna temperature, T* A , is reproduced for the 0-H2O 
1 10 — > loi line, observed by SWAS. However, the incorporated 
beam dilution factor is uncertain. The total signal is dominated 
by the signal of the PDR. 

Water line emission is found also outside of the molecular 
cloud core, but is a few orders of magnitude weaker than in 
the PDR, in agreement with observations (Snell et al., 2000). 
In this paper, freeze-out of water onto dust grains has not been 
taken into account. Bergin et al. (1995) probed the evolution of 
molecular clouds, including depletion of atoms and molecules 
onto grain surfaces in the temperature and density regime we 
are interested in. They find that 50-60% of the water vapour 
can be removed from the gas phase, when gas temperatures are 
20-25 K, dust temperatures 10 K and <nu> = 10 3 cm -3 , as in 
the case of the EMC. Consequently, depletion can result in an 
even lower intensity. Depletion of water molecules in the PDR 
will occur only inside the shielded clumps (not the irradiated 
clump edges). We find for nu ~ 10 5 cm -3 , T g = 25 K and = 
15 K that ~ 90% of the water vapour is frozen out (Bergin et 
al. 1995). 

Predictions have been made for the intensities and integrated 
antenna temperatures of various transitions that HIFI will probe 
in the future. It is found for S140 that the 2n — > loi and In — > 
Ooo transitions are the strongest for ortho- and para-H20, re- 
spectively, in all the models. 
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